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Abstract 



Method and apparatus for designing digital filters are provided. One method includes the 
steps of determining a real-valued discrete-frequency representation of a desired digital filter and 
transforming this discrete-frequency representation into a corresponding discrete-time representation. 
The discrete-time representation is circularly shifted and a window is applied to the discrete time 
representation to produce a zero-padded reduced length filter. Thereafter input signal x[n] is 
convolved with the reduced length filter. 



1 



DIGITAL FILTER DESIGN 

TECHNICAL FIELD 

The present invention relates to digital filtering, and in particular to digital 
filter design. 

BACKGROUND 

Some digital signal processing algorithms are based on filter design in the 
frequency domain. Typical applications are noise suppression algorithms for 
speech enhancement (see [1, 2]) and non-linear processors for echo cancel- 
lation. These filters are often quite long, which means that it may be desir- 
able to perform the convolution (filtering) in the frequency domain, since this 
is less complex than convolution in the time domain. Since the input signal 
to be filtered is typically much longer than the filter, the desired linear 
convolution has to be implemented by circular sectioned or block convolu- 
tion (see [3]) in order to avoid unacceptable delays and/ or complexity. 

A problem with filters designed in the frequency domain is that they are real- 
valued, which leads to a time domain representation in which the peak of the 
filter is split between the beginning and end of the filter (this is equivalent to 
a filter that is symmetric around lag 0, i.e. an non-causal filter). This makes 
the filter unsuitable for circular block convolution, since such a filter will 
generate temporal aliasing. This problem may be mitigated (but not elimi- 
nated) by weighting the data with a window that spans more than one block 
(see [2]). However, this will introduce a delay of 1 block, which is undesir- 
able. The split peak also makes the filter unsuitable for time domain convo- 
lution, since the important parts of the filter are at the beginning and end of 
the filter. This makes it difficult to approximate the filter with a shorter filter 
to reduce the complexity of the time domain convolution. 



Reference [4] describes a method that starts with an initial reduced length 
discrete frequency filter response (using either the Barlett or Welch method). 
The reduced length filter response is transformed to the time domain using a 
reduced length IDFT, circularly shifted to obtain a linear phase filter, zero- 
padded to obtain an extended length to avoid temporal aliasing, and trans- 
formed back to the frequency domain using an extended length DFT. A 
drawback of this method is that the resolution in frequency is reduced due 
to the initial reduced length filter. 

SUMMARY 

An object of the present invention is to provide a filter design method that 
avoids both temporal aliasing and the delay generated by multi-block 
windowing with improved frequency resolution. 

This object is achieved in accordance with the attached claims. 

Briefly, the present invention circularly shifts the filter in the time domain to 
recenter the filter peak and applies a short window around the peak to 
extract the essential part of the filter and automatically obtain zero-padding. 
This method has several advantages: 

1. The resulting zero-padded reduced length linear phase filter may be 
used to implement circular block convolution without temporal 
aliasing. 

2. If the window is sufficiently short, the convolution may also be per- 
formed in the time domain. 



3. If desired, the obtained linear phase filter may be transformed into a 
minimum phase filter to reduce the algorithmic processing delay. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



The invention, together with further objects and advantages thereof, may best 
be understood by making reference to the following description taken together 
5 with the accompanying drawings, in which: 

Fig. 1 is a time diagram illustrating a sample sequence x[n]; 

Fig. 2 is a time diagram illustrating a filter impulse response h[nj; 

Fig. 3 is a time diagram illustrating a time shifted filter impulse re- 
sponse h[n-l]; 

10 Fig. 4 is a time diagram illustrating a time shifted filter impulse re- 

sponse h[n-L]; 

Fig. 5 is a time diagram illustrating the convolution of sample sequence 
x[n] and filter impulse response h[n]; 

Fig. 6 is a time diagram illustrating a zero-padded sample sequence 

15 x[n]; 

Fig. 7 is a time diagram illustrating a zero-padded filter impulse re- 
sponse h[n]; 

Fig. 8 is a time diagram of a longer sample sequence x[n]; 

Fig. 9 is a time diagram of a first zero-padded block of sample sequence 
2 0 x[n] in fig. 8 padded with zeroes; 

Fig. 10 is a time diagram of a second zero-padded block of sample 
sequence x[n] in fig. 8; 

Fig. 1 1 is a time diagram of a third zero-padded block of sample se- 
quence x[n] in fig. 8; 

25 Fig. 12 is a time diagram of a diagram of the convolution between said 

first block in fig. 9 and the zero-padded filter impulse response in fig. 7; 

Fig. 13 is a time diagram of a diagram of the convolution between said 
second block in fig. 10 and the zero-padded filter impulse response in fig. 7; 

Fig. 14 is a time diagram of a diagram of the convolution between said 
30 third block in fig. 1 1 and the zero-padded filter impulse response in fig. 7; 

Fig. 15 is a frequency diagram illustrating a real-valued filter transfer 
function H[k]; 
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Fig. 16 is a time diagram of the filter impulse response h[n] that is 
obtained after transforming the transfer function in fig. 15 to the time domain; 

Fig. 17 is a time diagram of the filter impulse response of fig. 16 after a 
circular shift of N/2 samples; 
5 Fig. 18 is a time diagram illustrating a time window w[n]; 

Fig. 19 is a time diagram of the filter impulse response in fig. 17 after 
multiplication by the window in fig. 18; 

Fig. 20 is a time diagram of the filter impulse response of fig. 19 after a 
circular shift to remove leading zeroes; 
10 Fig. 21 is a frequency diagram of the magnitude of the filter transfer 

function F[k] that is obtained after transforming the filter impulse response of 
fig. 20 to the frequency domain; 

Fig. 22 is a frequency diagram of the phase of the filter transfer function 
F[k] that is obtained after transforming the filter impulse response of fig. 20 to 
1 5 the frequency domain; 

Fig. 23 is a time diagram of the filter impulse response of fig. 20 after 
transformation to a minimum phase filter; 

Fig. 24 is a flow chart of an exemplary embodiment of the convolution 
method in accordance with the present invention; and 
2 0 Fig. 25 is a block diagram of an exemplary embodiment of a convolution 

apparatus in accordance with the present invention. 

DETAILED DESCRIPTION 

25 Since the concepts linear and circular convolution are essential for the 

present invention, these concepts will be described in more detail with 
reference to fig. 1-14. 

Fig. 1 is a time diagram illustrating a short sample sequence x[n] of length L. 
30 This sequence is to be convolved with a filter having an impulse response h[n] 

illustrated in fig. 2. The definition of linear convolution in the time-domain is 
(assuming a linear time-invariant system): 



y[n] = x[ri] ® h[n] = ^x[k]h[n - k] = in this case = £ x[&]/z|> - k] 

k=-oo k=Q 



This expression may be interpreted in the following way: 

1. Form a set of shifted sequences {h[n], h[n-l], ... , h[n-(L-l)]} as illus- 
trated in fig. 2-4. 

2. Multiply each sequence by a corresponding sample x[0], x[l], x[L-l] 
of sequence x[n]. 

3. Add the scaled sequences to form the convolution as illustrated in fig. 5. 

From the above description and from fig. 5 it is clear that the resulting 
convolution y[n] will have a length of L+M-l samples. Thus, the filtered 
sequence y[n] will be longer than the original sequence x[n], which has only L 
samples. 

In the above description it was assumed that the convolution was performed 
in the time domain. However, it is also possible to perform convolution in the 
frequency domain by transforming each of the sequences x[n], h[n] to the 
frequency domain using the Discrete Fourier Transform, DFT (typically 
implemented by the Fast Fourier Transform, FFT), multiply the corresponding 
transforms with each other and perform an Inverse Discrete Fourier Trans- 
form, IDFT (typically implemented by the Inverse Fast Fourier Transform, 
IFFT), back to the time domain to obtain y[n]. Thus, in the frequency domain 
the convolution may be expressed as: 

Y[k] = X[k]H[k] 

This method is called circular convolution. In fact this method is often pre- 
ferred, since it reduces the complexity of the convolution, especially for longer 
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filters. However, there are some problems with this approach that have to be 
solved: 

1. Typically the input sequence x[n] and the filter sequence h[n] ore of 
5 different lengths, but the transformed sequences have to be of the same 

length to allow multiplication of the transforms. This problem may eas- 
ily be solved by zero-padding the shorter sequence up to the length of 
the longer sequence before the DFT operations. 

M 2. As noted above, the filtered sequence y[n] has a length of L+M-l sam- 

ples, while the input sequence x[n] has a length of only L samples. The 
first M-l samples of this output sequence are erroneous due to so called 
temporal aliasing caused by a wrap around of the last M-l samples in 
the output sequence (the last M-l samples in fig 5 would be added to 

15 the first M-l samples). The result is referred to as circular convolution 

with aliasing. However, what we want is a circular convolution that is 
equivalent to a linear convolution containing all L+M-l samples of y[n]. 
The proper way to handle this problem is to zero-pad both x[n] and h[n] 
up to at least the length N=L+M-1 of y[n] before performing the DFT, as 

2 0 illustrated in fig. 6 and 7. 

Fig. 8 is a time diagram of a longer sample sequence x[n]. The longer sequence 
may be filtered in the time domain using the definition of linear convolution in 
the time domain given above. Since the filter only influences a limited (M) 
25 input samples at a time, it is possible to calculate the output sample of y[n] 

"on the fly" as the filter h[n] is shifted through the input sequence x[n]. Thus, 
only the latest M samples of x[n] are required to calculate the current value of 
y[n]. However, if the filter is long the method is quite complex. In such cases a 
frequency domain method would be preferable. 

30 

In principle it would be possible to filter long sequences in the frequency 
domain in the way described above. However, practical considerations often 
make this approach unacceptable. For example, consider a sequence of 



7 



speech samples in a telephony application. Typically speech is sampled at a 
sampling rate of 8000 samples/ s. Thus, 1 second of speech will contain 8000 
samples. The complexity of the necessary DFT and IDFT would be consider- 
able for such large sequences. Another problem is that such an approach 
5 would lead to unacceptable delays, since the entire input sequence has to be 

collected before the DFT can be calculated. A better approach is to use block 
convolution in accordance with either the overlap-save or the overlap-add 
method. The overlap-add method will now be described with reference to fig. 
8-14. 

IX) 

In accordance with the overlap-add method the long input sequence x[n] of fig. 
8 is divided into blocks of length L as illustrated by fig. 9-11. Each block is 
zero-padded up to the required length L+M-l (where M is the filter length). 
Each zero-padded block is convolved with the zero-padded filter (se fig. 7) in 

15 the frequency domain as described with reference to fig 6-7.. This 

results in the convolved blocks illustrated in fig. 12-14. Finally, the 
convolved sequences are added to form the output sequence y[n]. It is 
noted that in the overlap regions indicated in fig. 12-14 the final result 
is obtained by adding partial results from two block convolutions. This is 

2 0 the reason for the name overlap- add. 

The overlap- save method is similar to the overlap-add method. However, 
in this case the blocks of the input sequence x[n] are overlapping (the 
last samples of one block are identical to the first samples of the next 
25 block). After frequency domain convolution, the erroneous first M-l 

samples of each filtered block (the samples that contain temporal 
aliazing) are discarded before the blocks are reassembled. 

In some applications the filter is determined in the frequency domain. For 
30 example, in telephony applications noise suppression based on spectral 

subtraction is often used (see [1, 2]). In this case the filter is determined as a 
function H(o) of frequency: 
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H(a>) 




where a, 0, 8 are constants and <f> v (» and 6, (a?) are estimates of the 



power spectral density of the pure noise and noisy speech, respectively. This 
expression is obtained from the model: 



where v[n] is the noise signal, x[n] is the noisy speech signal and y[n] is the 
desired speech signal. An estimate of the desired signal y[n], in which the 
noise has been suppressed, is obtained by applying the filter represented by 
H(o) to the noisy signal x[n]. 

Another example of an application in which the filter is determined in the 
frequency domain is a frequency selective non-linear processor for echo 
cancellation. In this case the filter is defined by the function: 



where <& x (a>) represent an estimate of the power spectral density of a signal 



spectral density of the echo signal e[n]. The filter is based on the model: 



x[n] = y[n] + v[n] 



x[n] contaminated by echo and O e (tf>) represents an estimate of the power 



x[n] = y[n] + e[n] 



where y[n] is the desired speech signal. An estimate of the desired signal y[n], 
in which the echo has been cancelled, is obtained by applying the filter 
represented by H(co) to the echo contaminated signal x[n]. 



In the above examples the filter is described by a real-valued continuous- 
frequency transfer function H(o). This function is sampled to obtain a discrete- 
frequency transfer function H[k]. This step is typically performed when the 
estimates are based on parametric estimation methods. However, it is also 
possible to obtain the discrete-frequency transfer function H[k] directly, for 
example by using periodogram based estimation methods. An advantage of 
parametric estimation methods is that the estimates typically have lower 
variance than estimates from periodogram based methods. 

Fig. 15 is a frequency diagram illustrating an example of a real- valued discrete- 
frequency filter transfer function H[k]. Fig. 16 is a time diagram of the filter 
impulse response h[n] that is obtained after transforming the transfer function 
in fig. 15 to the time domain using an IDFT. As may be seen from fig. 15 this 
filter has some unattractive features: 

1 . The filter has the full length N, which means that the already calculated 
transfer function H[k] of fig. 15 can not be used for frequency domain 
convolution without causing temporal aliasing as described above. 

2. If the convolution is performed in the time domain, the filter is not 
suitable, since its peak is split between the beginning and end of the 
filter h[n] . This will introduce a long algorithmic processing delay. 

In accordance with the present invention both these disadvantages may be 
avoided. This will now be explained with reference to fig. 17-22. 

The first step is to circularly shift the filter h[n] by N/2 samples. Mathemati- 
cally this may be expressed as: 



h[n] -> h[(n + N/2)modn] 



10 



Fig. 17 is a time diagram of the filter impulse response of fig. 16 after a 
circular shift of N/2 samples (N is assumed to be an even integer in this 
example). 

5 The next step is to apply a window to the shifted filter. Fig. 18 is a time 

diagram illustrating an example of a time window w[n]. In this case the 
window is a Kaiser window (see [5]) having a length of M samples (in the 
example M is assumed to be an odd integer less than N). Fig. 19 is a time 
diagram of the filter impulse response in fig. 17 after multiplication by the 
10 window in fig. 18. The result is a filter u[n] with leading and trailing zeroes. 

Preferably the leading zeroes are removed by another circular shift of D 
- : samples, where D is the number of leading zeroes. Fig. 20 is a time diagram of 

" = the filter impulse response of fig. 19 after a circular shift to remove leading 

:15 zeroes. The resulting filter f[n] will have its first non-zero tap at n=0 and will 

:; have N-M trailing zeroes. 

The result of the steps performed in fig. 17-20 is illustrated in fig. 21-22. Fig. 
2 1 is a frequency diagram of the magnitude of the filter transfer function F[k] 

2 0 that is obtained after transforming the filter impulse response f[n] of fig. 20 to 

the frequency domain. Fig. 22 is a corresponding frequency diagram of the 
phase of the filter transfer function F[k]. As may be seen from these figures the 
real-valued filter H[k] has been transformed into a linear phase filter F[k] 
(actually linear phase is obtained already after the first circular shift in fig. 

25 17). 

It is noted that the transformed filter f[n] in fig. 20 has an effective length of 
only M taps, and that the remaining N-M taps are zero. This fact may be 
exploited in two ways: 



1. The zeroes may be used to implement linear convolution in the fre- 
quency domain by using circular convolution of blocks of length L in 
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combination with the overlap-add method (fig. 8-14) or the overlap-save 
method. The block length L is obtained as: 

L<N-M + \ 

5 

2. The convolution may be performed in the time domain if the filter length 
M is selected sufficiently small. 

The GSM system will now be considered as an example of this procedure. This 
10 system already has a natural block length for block convolution, since speech 

is divided into frames of 160 samples. Thus, a natural choice is a convolution 
block length L of 160 samples. This leads to a natural FFT length N of 256 
samples. This results in a maximum filter length M of 97 samples (M<=N- 
L+l). 

W 5 

For a convolution block length of 256 samples the filter length typically lies in 
the interval 20-97 filter taps. For a Kaiser window the parameter p typically 
lies in the interval 2-5. Experiments have shown that a filter length M of 55 
samples obtained from a Kaiser window having a length of 55 samples and a 
2 0 parameter p«3 is a good choice, which even allows implementation of the 

convolution in the time domain using a signal processor. 

Fig. 20 illustrates that the linear phase filter f[n] is symmetric around the peak 
at (M-l)/2. This will lead to an algorithmic delay P of half the filter length. This 
25 algorithmic delay may be reduced by transforming the filter f[n] to a minimum 

phase filter. A minimum phase filter may be obtained from the linear phase 
filter in fig. 20 by using the method described in [6]. Briefly, this method starts 
with the expression: 



30 



/(*)= tfik]z k 

k=-n 
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where f[k] are the taps of a linear phase filter that is re-centered (linearly 
shifted) so that it is symmetric around tap 0 instead of tap (M-l)/2 as in fig. 
20. According to the spectral factorization lemma f(z) may also be written as: 

f(z) = Cg(z)g(z- 1 ) 

where C is a constant and g(z) is a polynomial 

g{z) = z n + g l z n ~ 1 +... + g n 

having all zeroes inside the unit circle. If g(z) is determined, a minimum 
phase filter with the same magnitude response in the frequency domain as 
the linear phase filter may be obtained as: 

Cg(z)g(V) 

Several methods exist for determining the coefficients of g(z) from the 
coefficients of f(z). For example, g(z) may be determined by solving the non- 
linear system of equations: 
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using the Newton-Rapson algorithm. Applying the spectral factorization 
lemma to the linear phase filter in fig. 20 results in a minimum phase filter 
with the same magnitude response (fig. 21) in the frequency domain. In the 
time domain the algorithmic delay P in fig. 22 has been reduced, as illus- 
trated in fig. 23. 
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In the above description of fig. 15-23 a certain processing order has been 
assumed. For example, the circular shift in fig. 17 was performed before the 
windowing in fig. 19. However, this order may be reversed if the window is 
circularly shifted accordingly. Furthermore, the circular shift in fig. 20 may 
be eliminated by including it in the circular shift in fig. 17 and circularly 
shifting the window accordingly. 

In this description a Kaiser window was used. However, all the common 
window types, such as Hamming, Hanning, Barlett windows etc. are 
feasible. Even a rectangular window may be used. 

Fig. 24 is a flow chart of an exemplary embodiment of the convolution method 
in accordance with the present invention. This embodiment relates to noise 
suppression and uses frequency domain convolution based on the overlap-add 
method and minimum phase filtering. Step SI collects the next input signal 
block. Using this block step S2 determines the transfer function H[k] from the 
above described noise suppression algorithm. Step S3 transforms H[k] to the 
time domain using the IFFT. Step S4 circularly shifts the time domain repre- 
sentation of the filter N/2 samples to produce a linear phase filter. Step S5 
applies a window to the linear phase filter to reduce its length. Step S6 
performs another circular shift to remove leading zeroes created by the 
window. Step S7 transforms the resulting reduced length linear phase filter 
into a minimum phase filter. Step S8 transforms the already zero-padded filter 
and the zero-padded input signal block to the frequency domain using the 
FFT. Step S9 multiplies the filter transform with the signal block transform to 
perform the convolution. Step S10 transforms the result back to the time 
domain using the IFFT. Step Sll combines the current convolved block with 
the previous block in accordance with fig. 12-14. Thereafter the algorithm 
returns to step SI to collect the next input signal block and repeat the 
procedure. 

If the convolution is performed in the time domain, steps S8- Sll are replaced 
by a time domain convolution step. 
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Fig. 25 is a block diagram of an exemplary embodiment of a convolution 
apparatus in accordance with the present invention. This embodiment is 
suitable for performing the method described in fig. 24. An input signal x[n] is 
forwarded to a buffer 10, which collects and stores a signal block. Using the 
input signal block a filter design block calculates the noise suppression filter 
H[k]. An IFFT block 14 transforms the filter to the time domain. A circular 
shift block 26 circularly shifts the time representation of the filter by N/2 
samples. A window block 18 applies a window to the shifted filter. Another 
circular shift block 20 removes leading zeroes from the reduced length filter. A 
minimum phase block 22 transforms the filter into a minimum phase filter. 
An FFT block transforms the filter back to the frequency domain. The trans- 
form is forwarded to a multiplication block 26. The input signal block from 
buffer 10 is also forwarded to a zero-pad block 20, which pads zeroes to the 
block up to the length of L+M- 1 . The zero-padded block is transformed to the 
frequency domain by an FFT block 30. The transform is forwarded to the other 
input of multiplication block 26. The convolved signal block from multiplica- 
tion block 30 is transformed back to the time domain in an IFFT block 32. The 
time representation of the convolved block is forwarded to an overlap buffer 34 
and to a combination block 36. Overlap buffer 34 extracts the overlap part 
(the last part) of the received signal block and outputs the stored overlap part 
of the previous convolved signal block. Combination block 36 adds the overlap 
part of the previous block to the beginning of the current block and outputs 
the remaining part of the current block unchanged. The result is the con- 
volved signal y[n]. 

Typically the blocks in fig. 25 are implemented by one or several microproc- 
essors or micro/ signal processor combinations. 

If the convolution is performed in the time domain, blocks 24- 36 are replaced 
by a time domain convolution block. 



Rather than repeating the above description, it is noted that a digital filter 
design apparatus in accordance with an exemplary embodiment of the 
present invention will include blocks 10-22 in fig. 25. 

It will be understood by those skilled in the art that various modifications 
and changes may be made to the present invention without departure from 
the scope thereof, which is defined by the appended claims. 
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U.S. CLAIMS 

1. A method of designing a digital filter, including the steps of 
determining a real-valued discrete-frequency representation of a desired 

full length digital filter; 

transforming said discrete-frequency representation into a correspond- 
ing discrete-time representation; 

circularly shifting said discrete-time representation; and 
applying a shortening window to said discrete-time representation to 
produce a zero-padded reduced length filter. 

2. The method of claim I, including the step of circularly shifting said 
reduced length filter to remove leading zeroes. 

3. The method of claim 1 or 2, wherein said real-valued discrete-frequency 
representation is formed by a noise suppressing spectral subtraction algo- 
rithm. 

4. The method of claim 1 or 2, wherein said real-valued discrete-frequency 
representation is formed by a frequency selective non-linear algorithm for echo 
cancellation. 

5. The method of claim 1, wherein said window is a Kaiser window. 

6. The method of claim 1, including the step of transforming said reduced 
length filter into a minimum phase filter. 

7. A digital convolution method, including the steps of 

determining a real-valued discrete-frequency representation of a desired 
full length digital filter; 

transforming said discrete-frequency representation into a correspond- 
ing discrete-time representation; 

circularly shifting said discrete-time representation; 
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applying a shortening window to said discrete-time representation to 
produce a zero-padded reduced length filter; and 

convolving an input signal with said zero-padded reduced length filter. 

8. The method of claim 7, including the step of circularly shifting said 
reduced length filter to remove leading zeroes. 

9. The method of claims 7, including the step of transforming said reduced 
length filter into a minimum phase filter. 

10. The method of claim 7, 8 or 9, including the step of performing the 
convolution step in the time domain using the discrete-time representation of 
said reduced length filter. 

11. The method of claim 7, 8 or 9, including the step of performing the 
convolution step in the frequency domain by using the overlap-add method. 

12. The method of claim 7, 8 or 9, including the step of performing the 
convolution step in the frequency domain by using the overlap-save method. 

13. A digital filter design apparatus, including 

means for determining a real-valued discrete-frequency representation 
of a desired full length digital filter; 

means for transforming said discrete-frequency representation into a 
corresponding discrete-time representation; 

means for circularly shifting said discrete-time representation; and 

means for applying a shortening window to said discrete-time repre- 
sentation to produce a zero-padded reduced length filter. 

14. The apparatus of claim 13, including means for circularly shifting said 
reduced length filter to remove leading zeroes. 
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15. The apparatus of claim 13 or 14, wherein said window applying means 
implements a Kaiser window. 

16. The apparatus of claim 13, including means for transforming said 
reduced length filter into a minimum phase filter. 

17. A digital convolution apparatus, including 

means for determining a real-valued discrete-frequency representation 
of a desired full length digital filter; 

means for transforming said discrete-frequency representation into a 
corresponding discrete-time representation; 

means for circularly shifting said discrete-time representation; 

means for applying a shortening window to said discrete-time repre- 
sentation to produce a zero-padded reduced length filter; and 

means for convolving an input signal with said zero-padded reduced 
length filter. 

18. The apparatus of claim 17, including means for circularly shifting said 
reduced length filter to remove leading zeroes. 

19. The apparatus of claims 17, including means for transforming said 
reduced length filter into a minimum phase filter. 

20. The apparatus of claim 17, 18 or 19, including means for performing 
the convolution step in the time domain using the discrete-time representation 
of said reduced length filter. 

21. The apparatus of claim 17, 18 or 19, including means for performing 
the convolution step in the frequency domain by using the overlap-add 
method. 

22. The method of claim 17, 18 or 19, including means for performing the 
convolution step in the frequency domain by using the overlap- save method. 
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ABSTRACT 

A digital convolution apparatus, includes means (12) for determining a real- 
valued discrete-frequency representation of a desired digital filter and means 
5 (14) for transforming said discrete-frequency representation into a corre- 

sponding discrete-time representation. Means (16) circularly shift the discrete- 
time representation and means (18) apply a window to the discrete-time 
representation to produce a zero-padded reduced length filter. Thereafter 
means (24-36) convolve an input signal x[n] with the reduced length filter. 
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